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We analyze and discuss the quantum noise in signal-recycled laser interferometer gravitational- 
wave detectors, such as Advanced LIGO, using a heterodyne readout scheme and taking into account 
the optomechanical dynamics. Contrary to homodyne detection, a heterodyne readout scheme can 
simultaneously measure more than one quadrature of the output field, providing an additional way 
of optimizing the interferometer sensitivity, but at the price of additional noise. Our analysis pro- 
vides the framework needed to evaluate whether a homodyne or heterodyne readout scheme is more 
optimal for second generation interferometers from an astrophysical point of view. As a more theo- 
retical outcome of our analysis, we show that as a consequence of the Heisenberg uncertainty princi- 
ple the heterodyne scheme cannot convert conventional interferometers into (broadband) quantum 
non-demolition interferometers. 

. PACS numbers: 04.80.Nn, 03.65.ta, 42.50.Dv, 95.55.Ym 
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I. INTRODUCTION 



Long-baseline laser-interferometer gravitational-wave (GW) detectors have begun operation in the United States 
(LIGO Q), Europe (VIRGO and GEO 600 |3]) and Japan (TAMA 300 H). Even as the first detectors begin the 
search for gravitational radiation, development of the next generation detectors, such as Advanced LIGO (or LIGO- 
II), is underway. With planned improvements in the seismic noise reduction — via active vibration isolation 
I and in the limits set by thermal noise — via the improved mechanical quality of the optics and clever suspension 
' strategies 0] , the sensitivity of second generation detectors is expected to be quantum- noise-limited in much of the 
, detection band from 10 to 10'' Hz. 

■ The optical configuration of all current GW detectors includes a Michelson interferometer. Two 4km-long Fabry- 
^ , Perot cavities are inserted into the arms of the Michelson interferometer; the optical field builds up in the cavities 
and samples the GW-induced phase shift multiple times. The arm cavities, thus, increase the phase sensitivity 
5^1^ _ of the detector. The Michelson-based optical configuration makes it natural to decompose the optical fields and 
bJ[)' the mechanical motion of the arm-cavity mirrors into modes that are either symmetric (i.e. equal amplitude) or 
^ ' antisymmetric (i.e. equal in magnitude but opposite in signs) in the two arms, as explained in detail, for example, in 
• 1— I , Refs. |3,HiS^3- No light leaves the interferometer from below the beam-splitter (BS) or dark port, except the lights 
induced by the antisymmetric motion of the arm-cavity test-mass mirrors, e.g., due to a passing-by gravitational wave, 
^ , or due to vacuum fluctuations that originally enter the interferometer from the dark port. Since GW interferometers 
- - ' operate on a dark fringe, the intensity of the light exiting the antisymmetric port is quadratic in the GW amplitude, 
and therefore insensitive to it, to first order. The standard way to circumvent this is to interfere the signal field with 
a relatively strong local oscillator (LO) field, such that the intensity of the total optical field, detected at the beat 
frequency, varies linearly with the GW amplitude. The various methods of measuring the GW-induced signal at the 
antisymmetric port are referred to as readout schemes. 

Previously 7. A. lol. Hol ITU, the quantum noise in Advanced LIGO was calculated assuming a homodyne readout 
scheme, in which the LO field oscillates at exactly the same frequency as the incident laser. The homodyne readout 
scheme can pose significant technical challenges for laser noise. In this paper we consider heterodyne readout schemes, 
in which the LO has different frequencies from the carrier. The heterodyne readout is usually implemented, as in 
initial LIGO (or LIGO-I), by using phase modulated light: the li ght incident on the interferometer consists of a 
carrier and radio frequency (RF) phase modulation (PM) sidebands |3]| . Using the Schnupp asymmetry 0| , the PM 
sidebands are transmitted to the photodetector as efficiently as possible, while the carrier still returns to the bright 
port. The transmitted sidebands then act as a LO against which the GW signal can beat. Demodulation at the 
modulation frequency converts the signal back down into the baseband. This technique circumvents laser technical 
noise by upconverting the signal detection to frequencies where the laser light is shot-noise-limited (a few MHz). 
Here we do not concern ourselves with technical noise on the laser; we consider only the fundamental quantum noise 
on the light. When the RF modulation-demodulation readout scheme is implemented, more than one quadrature of 
the interferometer output will be available for measurement, providing an additional tool for the optimization of the 
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FIG. 1: We draw a signal- (and power-) recycled LIGO interferometer. The laser light enters the interferometer from the left 
(bright port), through the power-recycling mirror (PRM), and gets split by a 50/50 beamsplitter into the two identical (in 
absence of gravitational waves) arm cavities. Each of the arm cavities is formed by the input test-mass (ITM) and the end 
test-mass (ETM) mirrors. A signal extracting mirror (SEM) is placed at the dark port, forming a signal extracting (SE) cavity 
with the ITMs. 



sensitivity, which is not available in homodyne detection. 

However, an additional quantum noise contribution, as compared with the homodyne readout scheme, usually 
appears in this scheme during the photodetection process — as was realized by Gea-Banacloche and Luechs in 
Ref. (12] where they evaluated the compatibility of squeezing and modulation-demodulation readout schemes in 
simple Michelson interferometers, and also by Schnupp, using more general considerations This additional 

contribution is due to vacuum fluctuations in frequency bands that are twice the modulation frequency away from 
the carrier. Subsequently, the heterodyne scheme was investigated in more detail by Niebauer et al. Jjj and Meers 
and Strain [l^ . These works 0i ITsL IT^ . ITsI focused exclusively on the detection of the output phase quadrature with 
phase modulated LO light (at the output port), which is appropriate for conventional GW interferometers with low 
circulating power, and hence negligible back action noise, but not for the advanced GW interferometers considered 
here. 

The main purpose of this paper is to further generahze the results obtained in Refs. 0,0,0,0], by including 
the possibility of detecting generic quadratures with LO light that are mixtures of phase and amplitude modulation 
to the carrier, and applying them to advanced GW interferometers, such as Advanced LIGO. In particular, we 
provide expressions and examples of the quantum noise, taking into account explicitly both the variable-quadrature 
optimization and the additional heterodyne noise. This lays the foundation for optimization of the detector sensitivity 
for specific astrophysical sources and for comparison between heterodyne and homodyne schemes from an astrophysical 
point of view. The results of these investigations will be reported elsewhere (IQ] . 

Recently, Somiya |0] showed independently the possibility of measuring different quadratures through heterodyne 
detection, and investigated the consequences for both conventional and signal-recycled interferometers. However, the 
additional heterodyne noise was not explicitly taken into account — it was hoped that, in certain sophisticated hetero- 
dyne schemes, the additional heterodyne noise becomes negligible, while the variable-quadrature optimization remains 
possible. However, as we show in this paper, the additional heterodyne noise is a direct consequence of the Heisenberg 
uncertainty principle, and will always exist as long as more than one quadrature is available for simultaneous measure- 
ment. Moreover, the Heisenberg uncertainty principle gives rise to a quantum limit to the additional heterodyne noise, 
which is frequency independent unless a frequency-dependent squeezing is implemented. This frequency-independent 
quantum limit will seriously constrain the power of the variable-quadrature optimization of heterodyne schemes in 
achieving (broadband) quantum non-demolition (QND) performances. In fact, for conventional interferometers, all 
quantum-limited heterodyne detection can be shown to be equivalent to a frequency-independent homodyne detection 
performed on an otherwise identical conventional interferometer with lower input laser power. 

This paper is organized as follows: in Section^ we describe the modulation/demodulation process and derive the 
demodulated output of signal-recycled interferometers in terms of quadrature operators and arbitrary heterodyne 
field amplitudes — taking the simplest sine- wave modulation/demodulation scheme as an example; in Section IlIII we 
derive the expressions of the quantum noise spectral density in this scheme, and apply them to initial and Advanced 
LIGO interferometers; in Section FlV Al we analyze a completely general modulation/demodulation scheme, derive 
a quantum limit for heterodyne measurements and discuss the consequence of this quantum limit for conventional 
interferometers. Finally, in Section we present our conclusions. 
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II. THE RADIO-FREQUENCY MODULATION-DEMODULATION SCHEME IN ADVANCED LIGO 

A. Overview of Advanced LIGO optical configuration 

The Michelson interferometer is operated on the dark fringe to minimize static laser power, and hence the shot 
noise associated with this hght, at the antisymmetric (dark) port. Since most of the Ught returns toward the laser, a 
partially transmitting mirror, the power- recyc ling mirror (PRM) is placed between the laser and the beam splitter to 
'recycle' the light back into the interferometer [l^ (see Fig.^. The optical configuration currently planned to achieve 
quantum- limited performance in Advanced LIGO uses the Resonant Sideband Extraction (RSE) technique in 
addition to power-recycling. In RSE, an additional partially transmitting mirror, the signal extraction mirror (SEM), 
is placed between the antisymmetric port of the beamsplitter and photodetector (see Fig^. 

The optical properties (reflectivity, loss) of this signal extraction mirror and its microscopic position (in fractions of 
the wavelength of the laser light, 1.064 /im) can significantly influence the frequency response of the interferometer 0, 
|20| . When the signal extraction cavity (SEC) — comprising the SEM and the input test-mass (ITM) mirrors of the 
arm cavities — is exactly resonant or anti-resonant at the laser frequency, the bandwidth of the entire detector can 
be increased or decreased by altering the reflectivity of the SEM. These two special cases are referred to as resonant 
sideband extraction (RSE) 19] and signal recychng (SR) ^20J, respectively. 

As the signal cavity is slightly offset (detuned) from resonance (RSE) or antiresonance (SR), the frequency at 
which the peak optical response of the detector occurs can be shifted to frequencies where other noise sources are 
not dominant. Note that, unlike conventional interferometers and tuned RSE/SR interferometers, the frequency 
responses of detuned configurations are no longer symmetric around the carrier frequency, with only one resonant 
peak located either higher or lower than the carrier frequency. As a consequence, although the interferometer will 
respond resonantly to GW's with a certain nonzero frequency, only one of the two (upper and lower) sidebands the 
GW generates symmetrically around the carrier frequency is on resonance. More generally, the upper and lower GW 
sidebands contribute asymmetrically to the total output field, which makes the GW signal appear simultaneously in 
both quadratures of the output field 0, ITol ITU . Detuned configurations are neither RSE nor SR in the original 
sense, but roughly speaking, such a configuration can be classified as either RSE or SR by looking at whether the 
bandwidth of the entire interferometer is broader or narrower than that of the arm cavity. Historically, since SR was 
invented earlier than RSE, some literature refers to all configurations with a signal mirror as "Signal Recycled". 

Since detuned RSE allows us to control the spectral response of the interferometer and optimize for specific as- 
trophysical sources, it has become a strongly favored candidate for Advanced LIGO jsj. A notable consequence is 
that with the high laser power of Advanced LIGO, the optomechanical coupling induced by detuned RSE/SR signifi- 
cantly modifies the dynamics of the interferometer, introducing an additional resonance at which the sensitivity also 
peaks 

B. Modulation and demodulation processes 

The RF modulation-demodulation scheme comprises two parts: the modulation-preparation and demodulation- 
readout processes. In this section we consider only the simplest case, sine wave modulation and demodulation. A 
more general discussion of modulation/demodulation schemes can be found in Sec. II VI 

Phase modulated light is incident on the interferometer. It is composed of the carrier at the laser frequency ujq ~ 
2 • 10^^ s~^, and a pair of phase modulation sidebands offset from the carrier by several MHz, so uj^^ ^ 2tt ■ 10^ s^^ ^ 
GW-sideband frequency < 27r • lO^'s"^. The detection port is kept as dark as possible for the carrier, while the PM 
sidebands at ojq ± Wm are coupled out as efficiently as possible to act as the local oscillator for the GW- induced carrier 
field that leaks out. Maximal RF sideband transmission is adjusted in two ways: (i) by a path difference in the arms of 
the Michelson that is arranged to be highly transmissive for the RF component of the field - the Schnupp asymmetry; 
and (ii) by matching the transmission of the power-recycling and signal-extracting mirrors so that the effective cavity 
comprising those two mirrors is critically coupled. In addition to the gravitational-wave readout, the PM sidebands 
are also useful for controlling the auxiliary degrees of freedom of the interferometer j2ll | . 

In Fig. 121 we show the outgoing optical field at the interferometer output in the frequency domain, which consists of 
the GW sidebands, the Schnupp sidebands and quantum fiuctuations at the output. The relative amplitudes of the 
RF sidebands are intentionally shown to be unequal. This is a case of unbalanced heterodyning, and is an intrinsic 
feature of the detuned RSE interferometer. As described above, when we move to detuned RSE, the SEC is detuned 
from perfect carrier resonance, such that the resonance peak of the signal cavity coincides with one of the GW signal 
sidebands, at the expense of the other one; consequently, the GW signal appears in both quadratures of the output 
field. At the same time, this phase shift in the signal cavity moves both RF sidebands off perfect resonance as well, 
which results in poor output coupling of both RF sidebands. This can be remedied by offsetting the RF sideband 
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FIG. 2: The outgoing field at the interferometer output consists of GW sideband signals around clio, the Schnupp sideband 
fields at UJQ ± uj-^a and quantum fiuctuations spread out at all frequencies. 

frequency - or conversely, the macroscopic length of the SEC - to make one of the RF sidebands resonant [2^ . 
Hence, detuned RSE leads to unbalanced heterodyne fields. Although the carrier is phase modulated before entering 
the interferometer, the heterodyne fields at the output port will no longer act as a pure phase modulation on the 
carrier. 

At the detection port, a standard heterodyne detection procedure is used to extract the GW signal: 



Outgoing light ^ Photodctection Mixing with cos(ijJni^ + 



Low — pass filter =^ Demodulated output 



The photodctection process consists of taking the square of the optical field shown in Fig. [3 This operation mixes the 
GW signal (and quantum fluctuations) located at frequency ~ wo with the RF sideband fields located at frequency 
Wo ± '^m- As a consequence the GW signal is measured in the RF band around ±Wm- By taking the product of (or 
mixing) the photodetection output with the demodulation function, cos(aJmt + 0d), the GW signal is down-converted 
back to low frequencies. The result is then filtered by a low-pass filter, yielding a frequency-independent quadrature 
that does depend on (pu. However, as we shall see more quantitatively in the following sections, in addition to the GW 
signal (and quantum fluctuations) centered at ~ ujq, quantum fluctuations at ^ luq± 2uj-a_-i also enter the demodulated 
output at the antisymmetric port. This gives rise to an additional noise term that is not present in a homodyne 
readout scheme. 



C. Demodulated output of LIGO interferometers 

The optical field coming out from the interferometer [see Fig. ^ can be written as a sum of two parts: 



The first term, 



L{t) = 



E{t) = L{t) + S{t) 



(1) 



(2) 



is the (classical) LO light composed of the Schnupp sideband fields at frequencies woi'^m, with (complex) amplitudes 
ZJ-i- and Z?_, respectively. The magnitude and phase of D± depend on the specific optical configuration. The two 
quadratures of the LO are generated by either amplitude modulation (first quadrature) or phase modulation (second 
quadrature) of the input light. The second term in Eq. J^l , 



S{t) 



dw 
27r 

+^ dn 



-A 
+A 



27r 

dn 
2^ 



[e-'"*5^+h.c.] , 
e-*('"°-'"»+'')*6^„_2^.„+ji + h.c. 
e-*('^°+^)*6^„+o+h.c. 
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(contributions at irrelevant frequency bands) , 



(3) 



contains both the (classical) GW signal and the quantum fluctuations of optical fields near ojq. Here A ^ Wm is the 
demodulation bandwidth |33j |. [For simplicity and clarity, we only write out explicitly the terms that will eventually 
contribute to the demodulated output.] 

The photocurrent from the photodetector is proportional to the square of the optical field, 

i{t)(xE^{t) = L^{t) + 2L{t)S{t) + S^{t), 

= [contributions at frequencies 0, ±2uji^, ±(2wo ± Swm)? from i^(t)] 
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[contributions at irrelevant frequency bands] 
[terms quadratic in fe, 6^ , from S'^(f)] . 



h.c. 



h.c. 



(4) 



After taking the product of (or mixing) i{t) with cos(ct;iu^ + </'d) smd applying a low-pass filter with cutoff frequency 
A, we obtain the demodulated output 



i{t) cos{ojmt + (Pb) '"°^=^'''''' 0{(j)D;t) 



f ^ [(i?+e-^- +7^_e'*°)e-'^*6.„+o + h.c.] 

J -A ^''^ 
J -A ^''^ 

^ [i?_e-*-e-"* &.„_2.„+n + h.c] , 

J -A ^''^ 



(5) 



where we assume the local oscillator to be strong enough that the quadratic terms in Eq. ^ can be ignored. It is 
convenient to recast the demodulated output (jsj in terms of quadrature operators by using the following relation (for 
A < Lu): 



-A 



dn 

2^ 



OJ + O 



dn 

2^ 



V2Aoe-'^'b'^+^in)+h.c. 



(6) 



Here A = Aq e*" {Aq, a G 5R) is an arbitrary complex amplitude, and the quadrature operator fo^_,_iL is defined as [see 
also Refs. |il3] 



where 



b-.in) 



6^ (r?) = 6^ (f7) sin C + ^2 (^) cos C : 



V2 



-b' 



V2i 



(7) 



(8) 



The superscript lo on the quadrature fields is added to emphasize that the quadratures are defined with respect to 
the central frequency to. By applying relation © to the demodulated output lO), we get 



0(</)D;i) - 

in which we have defined 



+A 



dn 

2^' 



-int 



V2 



Do b'^^in) + \D+\ 6^«+2"-(17) + ILI-I 6^»-^'^»(J7) 



h.c. 



(9) 



(10) 
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and 

Co = I + arg {D+ e-^^- + e^^°) , (11) 
C± = ±</)D + | + argi5±. (12) 

In the frequency domain, we have, 
0((/)D;r!) = V2Do 

The first term inside the parenthesis, 6"° is an output quadrature field around the carrier frequency wq, which 
contains both the GW signal and vacuum fluctuations in the optical fields near the carrier frequency. In Refs. P, 
1^, llCt ITU , this quadrature field is related to the input quadrature field at the antisymmetric port via the input- 
output relations, from which the spectral density of the quantum noise can be derived. Measuring this field is 
the task of all readout schemes. For example, a homodyne scheme can measure directly an arbitrary frequency- 
independent quadrature. For this reason, we call the quadrature field 6^° the homodyne quadrature for distinction. 
The two additional terms inside the parenthesis are the additional noise, which come from vacuum fiuctuations around 
luq ± 2iOm- The sum of all three terms is what we measure in the heterodyne scheme, which we call the heterodyne 
quadrature. 



Quantity 


Symbol and Value 


Laser frequency 


ujo = 1.8 X lO^'' sec"^ 


GW sideband frequency 


n 


Input test-mass transmissivity 


T 


Arm-cavity length 


L = 4km 


Mirror mass 


m 


Light power at beamsplitter 


lo 


SEM amplitude reflectivity and transmissivity 




SEC length 


/ ~ 10 m 


SEC detuning 


(j) = [l^oVc]mod2ir 



TABLE L Basic quantities of Advanced LIGO interferometers 



1^ 

Dn 



Dn 



m < A, 



(13) 



D. Features of the RF modulation-demodulation scheme 



As can be inferred from Eq. (|ll|l . as long as \D+\ ^ \D-\, all homodyne quadratures can be measured through 
some heterodyne quadrature with the appropriate demodulation phase 0d- The (single-sided) spectral density S{fl) 
associated with the noise /i„ can be computed by the formula [see Eq. (22) of Ref. Q]: 

2TrS{n - n') S{n) = (in|/i„(r2) hl{fl') + hl{n') hn{n)\m) , (14) 

and if the input state of the whole interferometer is the vacuum state (|in) = |0a)), the following relation holds: 

{Oa\a^{n)al{n') + al{n')a,{n)\Oa) = 2TTS{n-n')s,j. (15) 

From Eq. H13() we see that the noise spectral density in the heterodyne quadrature is a sum of that of the homodyne 
quadrature, S]^°"'{(j)B-M), and those of the additional noise terms, ^^^'^^(Vd; f^)- Since 5^'°'°((/)d; f^) and S^'^'^{<j)B-M) 
come from different frequency bands, we assume that they arc uncorrelated, hence 

5^'(0d; ^^) = ^^■"(<^d; ^^) + St'^^i^j.^n) . (16) 

Assuming that the fields associated with the additional heterodyne noise are in the vacuum state, we get a white 
(frequency-independent) spectrum for the additional noise. 



\D+\ ^ + \D-\' ^ \D+\' + \D-\' 



(0d) - f^a - in „-,.d,r. , n „id,r.\1 ' (1 
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which usually depends on (j^D, unless either D_|_ or D_ is zero, which we refer to as the totally unbalanced case. In the 
case of balanced modulation, when \D^\ — |-D_|, only one quadrature, 



is measured, with additional noise 



^balanced = | + i (arg D+ + arg D_) , (18) 

^add balanced _J_ (19) 

and with a frequency-independent optimal demodulation phase 

0^^'^"'="'* = i[argi:i+-argD_]+iV7r, = 0, ±1, . . . . (20) 

This is the lowest possible additional noise for heterodyne schemes with just one pair of sidebands. The noise spectral 
density can have different shapes as a function of the homodyne angle At different signal sideband frequencies, 

the optimal homodyne angle Copt that gives the lowest homodyne noise can be different. In homodyne detection, 
since both quadratures of the carrier are generally not available, only a single frequency-independent quadrature can 
be measured (s^. By contrast, in heterodyne detection schemes (except for the balanced case), all quadratures are 
available for simultaneous measurement, and the final heterodyne noise at each frequency will be the minimum of all 
quadratures. 

III. NOISE SPECTRAL DENSITY AND THE EFFECT OF THE ADDITIONAL NOISE 

In this section, we write down the noise spectral density for both conventional and RSE interferometers when the 
RF modulation-demodulation scheme described in Section Mis used. 



Symbol 


Quantity 


Expression 


7 


Half bandwidth of arm cavity 


Tc 
4L 


13 


Phase gained by resonant field in arm cavity 


arctan(f2/7) 


hsQL 


Free-mass standard quantum limit 


J 

V mn2L2 


/sQL 


Characteristic input power for conventional 
interferometer to reach the SQL at f2 = 7 




K. 


Coupling constant for radiation-pressure effects 


2(/o//sql)7' 

^2(^2 +^2) 



TABLE II: Quantities derived from those listed in Tab.Q 



A. Total noise spectral density 

The input-output relation for RSE interferometers, including optomechanical effects, were derived in Refs. 0, |^ 
[see Eqs. (2.20)-(2.24) and (2.26) of Ref. ,^J]. The output fields in the frequency band of {luq — A, ujq + A) are (in the 
conventions used in this manuscript) 



6^° \ ^ 
62 ° / M 



C21 C22 / V "^2 / V -^2 y ft-sQL 



(21) 



where 



M = l + p^ e^''^ - 2p e^''^ ( cos 20 + ^ sin 20 ) , (22) 




FIG. 3: In the left panel we show the square root of noise spectral density, in units of hsQhil), for a conventional interferometer 
with 7o = /sQL, using balanced heterodyne modulation (solid curve) and the homodyne (dashed curve) scheme, plotted as 
functions of Q/^. The second quadrature is measured. In the right panel we plot the noise spectral density of the same 
interferometer, using unbalanced heterodyne detection, with homodyne angle chosen at the optimal value for the homodyne 
case, ({il) = Copthom(f2) (dashed line), and at the re-optimized value for the heterodyne case ~ Copthot(^^) (solid curve), 
respectively. The optimal heterodyne noise spectral density without the additional noise is also shown (dotted curve), which 
agrees with the result for frequency-independent homodyne detection Q. The SQL line is shown in both panels as gray straight 
lines. 

Cii ^ C22 = {l + p^) (^cos 20 + ^ sin 20^ - 2p cos 2(3 , (23) 
C12 = -T^ (sin 2(1) + IC sin^ 0) , C21 = (sin 2^ - JC cos^ 0) , (24) 



Di = + sin 



-i-l + pe^'^) cos 



(25) 



The quantities /C, /3, p, 0, r and /isql are defined in the same way as in Refs. 8, 0, 113 • We denote by h{VL) the 
gravitational strain and give a summary of the main quantities in Tables |2 and We assume that the fields a'^° 
incident on the unused input of the antisymmetric port are in the vacuum state for all frequencies. Moreover, the 
additional heterodyne noise fields J)'^o±2w„ (jl^ji must also be in vacuum states, since they are far away from the 

carrier frequency and are not affected by the ponderomotive squeezing effects of the interferometer. We assume that 
the higher-order terms of the modulation are not resonant in the interferometer, which is in general the case. Even 
if the higher-order sidebands are resonance, we would not expect any ponderomotive squeezing since the frequency 
is too high for the test-mass displacement to respond to an external force. Using Eqs. I|14|) - (|17|l and (|21|l . we obtain 
the total heterodyne noise spectral density in h, as a sum of the corresponding homodyne noise (first term) and the 
additional heterodyne noise (second term) [see Eqs. H1U|I . for the definitions of Z?o and ^0]: 



/)2 



"■SQL 



1 



2/C r2 \Di sin Co + D2 cos Co I 
1 



2/C r2 \Di sin Co + D2 cos Co I 



(Cii sin Co + C21 cosCo)^ + (C12 sin Co + C22 cosCo)^ 



\D+\' + \D\'' 



\D+ e-'^-o + £)_ e^^Dl 



(26) 



We note that the optimal heterodyne noise spectral density at a given GW signal sideband frequency is the minimum 
of those obtained varying (/)d (and thus Co)- 



B. Conventional interferometers 



For the power-recycled Fabry-Perot Michelson optical configuration, the so-called conventional interferometer, the 
GW signal appears only in the second (or phase) quadrature. Furthermore, barring imperfections, the transmission 
of the Schnupp sidebands is balanced. In our notation such a scheme is obtained by setting £)_ = —D*^, with 



9 



(pD = ±7r/2 + argD_|_ which is the optimal demodulation phase for all frequencies [see Eq. (|20ll ]. Evaluating Eq. H26|) 
in the case cj) — 0,t — 1, we get 

5r™"" = %^(^'+i+^) , (27) 

where the last term inside the parenthesis is the additional heterodyne noise, which is equal to 1/2 the shot noise 
in homodyne readout scheme (second term), originally derived in Ref. 15j. In the left panel of Fig. wc plot the 
noise curves of a conventional interferometer with Jq = /sql, using homodyne and balanced heterodyne detection, 
respectively, with the second quadrature measured. This is exactly the result in Ref. 0|. More sophisticated 
modulation schemes that can further lower or eliminate the additional heterodyne noise in this quadrature have been 
investigated by Schnupp 13], Niebauer et al. [iJl, and Meers and Strain 15]. 

If, on the contrary, the RF sidebands at the antisymmetric port are not balanced, one can measure arbitrary quadra- 
tures by adjusting the demodulation phase [see Sec. Ill P) . As proposed by Vyatchanin, Matsko and Zubova [23], and 
further investigated by Kimble, Levin, Matsko, Thorne and Vyatchanin (KLMTV) 7], measuring different quadratures 
at different GW signal sideband frequencies can allow conventional interferometers to beat the Standard Quantum 
Limit (SQL) [23] significantly, thus converting them into a QND interferometer. Somiya proposed that, by using 
a frequency-dependent demodulation phase, a KLMTV-type, frequency-dependent optimization is achievable in a to- 
tally unbalanced modulation scheme. However, the effect of the additional heterodyne noise was not explicitly taken 
into account and we show in this section that the additional heterodyne noise plays an important role as soon as one 
approaches the SQL. So much so, that for totally unbalanced heterodyne detection, the SQL cannot be beaten, and 
for intermediate levels of imbalance the SQL is beaten by very modest amounts. 

For simplicity, we first consider a totally unbalanced modulation scheme (which was the case investigated by 
Somiya ^17J). in which only (or only D_) is non-zero. From Eq. (|26|l . fixing r = 1, p = and <^ = 0, we have 



/j2 

ohct conv _ _ "-SQL 



(/C-tanCo)^ + l+ ^ 



cos2 Co 



(28) 



where the last term inside the parenthesis is the additional noise due to heterodyne detection. Using the optimal 
detection angle in the (frequency dependent) homodyne case QiE^, Copthom = arctan/C, one has 

Sh='^{lC'+2)>V2hl^^, (29) 
which cannot reach the SQL. Re-optimizing the detection angle, we obtain Copthet = arctan(/C/2). This gives 

[^r™"1.p, = %(^ + 2)>/^iQ,, (30) 

which only touches, but never beats, the SQL. In the right panel of Fig.|31 we plot the noise curve of a conventional 
interferometer with Iq = IsQh, the heterodyne noise spectral density using Copthom [given by Eq. (|29|l ]. and the 
optimal heterodyne noise spectral density [given by Eq. Ij^Uf) ]. As can be further verified, having two sidebands with 
unequal amplitude can allow the interferometer to beat the SQL, but only by very moderate amounts, and in limited 
frequency bands. 

We might still expect to use more sophisticated modulation-demodulation schemes to lower the additional hetero- 
dyne noise while retaining the possibility of variable-quadrature optimization. However, as we shall see in Sec. IIVI 
such an effort will be significantly limited by the Heisenberg uncertainty principle. 



C. Signal-recycled interferometers 



In this section, we give some examples of noise curves of detuned RSE interferometers with a heterodyne readout 
scheme, and compare them to the homodyne cases. 

In the balanced scheme the additional heterodyne noise is the lowest, but only one quadrature can be measured. 
For this case we show the effect of the additional heterodyne noise on the sensitivity curves in Fig.^ In the left panel, 
we plot the noise curves for a detuned RSE interferometer with T = 0.033, p = 0.9, 4> = — 0.47, Iq = /sql and 
m = 30kg, (the configuration considered in Refs. lESEl) when the first (C = 7r/2) and second quadratures (C — 0) 
are measured, by homodyne and balanced heterodyne read-out schemes. In the right panel, we plot the ratio of the 
heterodyne noise to the corresponding homodyne ones. The additional heterodyne noise has more features around 




FIG. 4: A detuned RSE interferometer (T = 0.033, p = 0.9, (j) = 7r/2 - 0.47, m = 30 kg, 7o = /sql) using balanced heterodyne 
detection. In the left panel we plot the square root of total heterodyne noise spectral density in the first (dash-dot curve) and 
second (long dashed curve) quadratures, compared with the homodyne ones (dot curve and solid curve, respectively) . The SQL 
line is also shown as a gray straight line. In the right panel we show the ratio of the square roots of the heterodyne and the 
homodyne noise spectral densities, for the first (dashed curve) and second (solid curve) quadratures. 



the two valleys of the noise curves, where the optomechanical dynamics (the RSE transfer function) determines the 
shape of the curves. Above ~ 200 Hz, the ratio between the square roots of the heterodyne and the homodyne noise 
spectral densities assumes the constant value ■y/3/2 ~ 1.22, which is due to the additional heterodyne noise when the 
shot noise dominates [see Eq. (I?7|l ]. 

Practical implementation of the RF sidebands in the interferometer has shown that detuned RSE configurations 
are likely to be very unbalanced 0, |22| ■ In the left panel of Fig. O we plot the unbalanced heterodyne noise spectral 
densities for the same interferometer parameters used in Fig. ^ with Co = 0, 7r/4, 7r/2 and 37r/4, and the optimal 
heterodyne noise obtained by maximizing over Co at each sideband frequency. Indeed, in the heterodyne readout 
scheme we have the advantage of optimizing the detection angle at different frequencies. At each particular signal 
sideband frequency, the optimal heterodyne noise spectral density is just the minimum of all quadratures. In the 
right panel of Fig. [Sj we compare the optimal heterodyne noise with the homodyne noise at C = and C, = 7r/2. 
As we see from this example, for the same interferometer configuration, neither the homodyne nor the heterodyne 
readout can provide a noise spectral density that is the lowest for all GW signal sideband frequencies. To make a 
more rigorous comparison between these two schemes a more critical study is required that takes into consideration 
specific astrophysical GW sources, the experimental feasibility and the other sources of noise, as well. [As an example, 
current Advanced LIGO design estimates the dominant, thermoelastic component at about the SQL |2J|. To lower the 
thermoelastic contribution below the SQL an interesting and challenging proposal has been analyzed recently p^.] 
The optimization of homodyne versus heterodyne readout schemes which include those effects is currently underway, 
and will be reported elsewhere [l6| . 

As shown in Ref. 0|, detuned RSE interferometers have an unstable optomechanical resonance. In the parameter 
regime emphasized in Refs. [ili[l3|ll|, the unstable resonance lies within the observation band — which gives a dip 
in the noise spectrum. Consequently, the control scheme must sense and act on the motion of the system within the 
observation band. In Ref. 0, an idealized control scheme is conceived for the homodyne readout, which suppresses 
the instability and leaves the noise spectral density unchanged. The same control issue will need to be addressed with 
the heterodyne readout scheme as well. 



IV. MORE GENERAL DISCUSSION OF HETERODYNE SCHEMES: MINIMAL ADDITIONAL NOISE 

AND QUANTUM LIMIT 

In Sec.|n]we discussed the sinusoidal modulation-demodulation scheme, which is the easiest to im plem ent. There 
exist more sophisticated schemes, such as those proposed by Schnupp, investigated by Niebauer et al. [lj|, and Meers 
and Strain ^15i |. that can further optimize the interferometer performances. These authors restricted their analyses 
to low-power interferometers and focused on the detection of the second (or phase) quadrature. In this section, we 
extend their discussions to the more general case where all quadratures can be measured. As we shall see, although 
modulation/demodulation readout schemes offers the advantage of variable-quadrature optimization, they are in 




FIG. 5: A detuned RSE interferometer (T = 0.033, p = 0.9, = tt/2 - 0.47, m = 30kg, Jo = Isql, same as Fig. using 
totally unbalanced heterodyne detection. In the left panel we show the noise curves for quadratures with C = (solid curve), 
7r/4 (long dashed curve), n/2 (short dashed curve) and 37r/4 (dash-dot curve) are shown, together with the final heterodyne 
noise optimized at all sideband frequencies (thick solid curve). In the right panel, the final heterodyne noise (thick solid curve) 
is shown along with the homodyne noise in the first (dashed curve) and second (dash-dot curve) quadratures. The optimal 
heterodyne noise without the additional noise is also shown (dotted curve) for comparison. The SQL is plotted in both panels 
as gray straight lines. 



general limited in converting non-QND interferometers to (broadband) QND interferometers. 



A. Quantum Limit for the additional heterodyne noise 



The field coming out from the dark port can be written, in time domain, as 

E{t) = [A{t) cos ojnt + P{t) sin uJot] + [Ei{t) cosojnt + E2{t) sinuJot] , (31) 

where the first term is the transmitted Schnupp sideband fields in the form of a combination of amplitude modulation 
[A{t)] and phase modulation [P{t)] to the carrier. In Eq. H31|l we denoted by Ei{t) and E2{t) the quadrature fields 
containing GW signal and quantum fluctuations. The output from the photodetector is then 

t{t) oc A{t)Ei{t) + P{t)E2{t) . (32) 

The amplitude and phase modulation is in general periodic functions, with the same angular frequency uj^ [wo 3> 

A{t) = ^Afce-^*^""*, Ak^AU: (33) 
fe 

P(t) - Pke-'"^--' , Pk = Plk ■ (34) 
fe 

(35) 

In the frequency domain Eq. H32|l reads: 

cx [Afefc^" {n - kujn,) + PkhT - fct^m)] • (36) 

k 

Denoting the demodulation function with D{t), the demodulated output is: 

0{t) = D{t)i{t) cx D{t)A{t)Ei{t) + D{t)P{t)E2{t) . (37) 
The demodulation function D{t) should have the same frequency as the modulation functions, therefore: 

D{t) = YDke-'''^"^\ Dk=D*_^. (38) 

k 
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Note that Eq. (|37|) is a generalization of Eq. (4) of Ref. Using the above equations, the Fourier transform of the 
demodulated output (|37|l can be written as: 

d{n) = [AkbT m + Pkb'^"m +Y^Y.Dl [^fe+p^r (^^ - P^™) + Pk+pbT {n - pu;™)] . (39) 

k p^O k 

Let us suppose that the low-frequency component of 0{il) is filtered out, then the first term in Eq. H39(l gives a 
frequency- independent quadrature field near wg, while the second term gives the additional heterodyne noise that 
arises from quantum fluctuations near loq ± puj^, with p = ±1, ±2, ±3, .... Since Wm ^ ^^gw, these fields are 
not affected by ponderomotive squeezing effects in the interferometer arm cavities and will be in the vacuum state. 
As a consequence, the additional heterodyne noise will also be frequency independent. [Unless frequency-dependent 
squeezed states are injected into the dark port of the interferometer.] In this way, for any particular quadrature C, 
there is a uniform minimum of the additional heterodyne noise at all frequencies. 

Let us now construct for an arbitrary quadrature C the optimal demodulation function D{t) and evaluate the 
minimal additional noise. 



or in the time domain. 



If we want to measure 6^° , Eq. (|39(l says that we have to impose 
lY.DlAk,Y.DlpA -(sinCcosC), 

\ fc k ) 



(sin C, cos C)=(^^£ D{t)A{t)dt, 1 £ D{t)P{t)dtj 



(40) 



(41) 



where T = 2Tr/L0m is the common period of the modulation and demodulation functions. Note that, in order for the 
quadrature ( to be measured, Eqs. l|in|l and l|lT|l need only be true up to a constant factor. Having written them in 
the current way, we have in fact chosen a specific normalization for D{t). Using the Parseval theorem and Eq. 141|1 . 
we derive for the spectral density of the additional heterodyne noise: 



^add _ ^ 




2 

+E 


Y.^kPk+p 


2 






p#o 

2 


k 


2 


= E 


E ^*kAk+p 

k 


+E 

p 


E ^*kPk+p 
k 





Y.^kAk 



Y.^iPk 



^ j\\t)[A^{t) + P\t)\dt- j\{t)A{t)d^ j\{t)P{t)dt 
^[ D\t) [A\t) + P\t)\dt-l. 



(42) 



[Note that Eq. (|^ . which corresponds directly to Eq. (12) of Ref. |0, is also consistent with Eq. (18) of Ref. [T5l| . 
since Eqs. H4()|l and (|41(l have already imposed a normalization for D{t)] In order to find the D{t) that satisfies 
Eq. H41|) and minimize S^'^'^ , we introduce two Lagrange multipliers, A and /i, and impose, 



5 J dt^ [{A^it) + P'^it)] D^{t) - 2XA{t)D{t) - 2tiP{t)D{t)'j =0, (43) 

which yields: 

XA{t)+^,Pit) 

^' A^t) + P^t)- ^ ' 

[In Eq. H43|l . the factors of 2 in front of A and ^ are added for simplicity.] Inserting Eq. (|44|l back into Eq. 141|l gives, 

A \ / sinC 



M( J = icosC )' ^''^ 
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where 



M = 



T 



T.,0 



\tJ, A^t) + p^t) 



A^{t) + P^{t) 

PHt) 



dt 



(46) 



A^t)+P^t) I 



The optimal demodulation function for the C quadrature is then given by inverting Eq. (|45|l and inserting the resulting 
A and into Eq. (|44ll . The minimal additional noise can be then obtained by inserting the optimal demodulation 
function into Eq. Wh : 



^addmin(^) = ( A A* ) M ^ ) - 1 = ( siu C COS C ) ^ ^'^"^ ^ ) - 1 = ( siu C COS C ) {M.'^ I) ^^^^ ) . (47) 

Moreover, we note an interesting property of M: 

I- M = (detM)M-^ (48) 

As a consequence, 

- I = M-i(I -M) = (detM)(M-i)2, (49) 



so 



det [M-i -I] = 1. (50) 

This implies that, the minimal additional noise can be written in a form 

oaddmin.^N ^ „• ^ N ( COS sin \ ( \( COS0 -sin0\ ( sinCA 

b (C) - ( smC cose ) I ^^^^ II e-fl ) 1 gi^^ cos0 ) \ cosC j ' ^^^^ 



with and i? frequency- independent, and determined by the eigenvectors and eigenvalues of the matrix M^^ — 1, 
which are determined ultimately by the amplitude and phase modulations. It is interesting to note that this minimal 
noise spectrum is of exactly the same form as that of a squeezed state. 

This phenomenon could in fact be anticipated from quantum mechanics. For the same sideband frequency f2, the 
different quadratures do not commute with each other, and have the following commutation relations: 

= 27ri sin(C - O^iS^ - ^) ■ (52) 

As a consequence, quantum fluctuations in the various quadratures are constrained by the Heisenberg uncertainty 
principle. As is well known, the squeezed states have the minimum noise spectrum allowed by the uncertainty principle. 
In modulation/demodulation schemes, all quadratures can be read out, with additional noise: 

6'^<=*(C,r!) =&^''(r!)+n(C,rJ). (53) 

So, all output observables should commute with each other, and as a consequence 

[6'^^*(C,ri),6^"*^(C',f^')] =0. (54) 

Since 6"° ($7) and ri(C, $7) come from different frequency bands of the output field, they must commute with each other, 
so we must have that the mutual commutators of n(C, fi) cancel those of fe^''(il): 

[n(C,f7),nt(C',17')] =-27risin(C-C')'^(^^-^^')- (55) 

Since they do not commute with each other, the additional noise n(C, £7) is also subject to the constraint of the 
Heisenberg uncertainty principle — in the same way as 6^°(f2), since the commutators only differ by a sign [see 
Eqs. (|52|l and H55|l ]. This explains why the minimum additional heterodyne noise has a spectral density of the 
same form as the squeezed states. The minimum noise spectrum (|51|) can be regarded as a quantum limit for 
modulation/demodulation schemes . 



&^°(17),6^°t(f7') 
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B. Impact of the Quantum Limit on conventional interferometers 



As discussed in Refs. |7|,|23||, using an appropriate readout scheme, conventional interferometers can achieve QND 
performance through a cancellation between shot and radiation-pressure noises. If the quadrature ^ is measured, we 
have 



Lhet conv 



(C,l^) =cosC 



e^'^ (tanC-/C)a^° + 



+ e 



i/3 



2/C- 



^SQL 



(56) 



and if we choose to measure the quadrature with C = Cc = arctan/C [35l | the part of the shot noise [the term proportional 
to tan C inside the bracket of Eq. (|56|l ] cancels the radiation-pressure noise [the term proportional to /Ca^" inside 
the bracket of Eq. (|56|l ] . The remaining shot noise [obtained from the term proportional to 0,2° inside the bracket 
of Eq. I|56fl ]. normalized to unit signal strength, is inversely proportional to /C, and it can be made lower (eventually 
lower than the SQL noise) by taking larger /C. However, for larger values of /C, tan<^c grows and the corresponding 
cos^c decreases. As can be seen from Eq. (|56|l . this implies an even smaller signal strength in the detected quadrature, 
which makes the additional noise, n(Cc,^^), more and more important. In fact, more generally the additional noise 
limits the extent to which the interferometer can beat the SQL. Writing the total heterodyne noise spectral density 
[of which Eq. (|28|) is a special case] , as 



Sh — 



"■SQL 

2JC 



(/C-tanC)^ + l + 



cos^ C 



(57) 



and following the argument that had led us to Eq. (|30|l . we obtain the following lower limit for the heterodyne noise: 



(58) 



where Copt(f^) is the optimal detection quadrature at frequency H., which depends also on the shape of S^'^'^{Q. Equa- 
tion H58|l says that, in order to beat the SQL significantly, the additional heterodyne noise at the optimal quadrature 
has to be much smaller than unity. However, since the additional heterodyne noise is frequency independent, and 
subject to the quantum limit (|51|l . this requirement cannot always be fulfilled if the optimal homodyne quadrature 
varies significantly with frequency in the observation band. As a consequence, heterodyne schemes will have very 
limited power in converting conventional interferometers into (broadband) QND interferometers. 

Due to the simplicity of the input-output relations of conventional interferometers, we can go a step further and 
obtain a cleaner result in this case. Let us suppose that the additional heterodyne noise have exactly the form of 
Eq. H51II . with generic values of and e^, i.e. it is quantum limited. Inserting Eq. (|51|l into Eq. H57II . we find the 
frequency-dependent optimal detection phase, 



tan Copt (f^) 



[1 - tanh2(i?/2)] IC{n) + 2 tanh(i?/2) sin 20 
2-h2tanh(i?/2)cos20 



and obtain: 

^quantum = 



[1 + 2cos(20) tanh(i?/2) + tanh2(i?/2)]/C2(r2) -4sin(2<?!)) tanh(i?/2) A:(f7) 
4/C(r2) [1 -I- cos(20) tanh(i?/2)] 



(59) 



(60) 



Moreover, the quantum- limited heterodyne noise spectral density (|60|l can be recast into exactly the same form as 
that of a frequency-independent homodyne detection 



^quantum = 



2/Ccff(f^) 



[(/Ceff(r!)-tanCcff)' + l] , 



(61) 



with 



tan Ccff 



tanh(i?/2) sin(2(/)) 
1 + tanh(i?y2) cos(20) 



/Ceff(51) = 



1 - 



l-tanh^(i?/2) 
2tanh(i?/2)cos(20) 



(62) 



Note that in the definition of JC^b the quantity multiplying IC (which is less than 1, since —1 < tanh(i?/2) < +1) 
can be absorbed into the input power [see the definition of IC in Table Equations (|61|l and (|62|l therefore 
relate a conventional interferometer with a quantum limited heterodyne readout scheme to an identical conventional 
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interferometer, but with lower input power and a frequency-independent homodyne readout scheme. As discussed by 
KLMTV, the latter does not exhibit broadband QND behaviour [although fine tunings of parameters can sometimes 
give a moderate SQL-beating noise spectral density]. This means, that the variable-quadrature optimization provided 
by heterodyne readout schemes does not enhance the QND performance of conventional interferometers at all. 

Nevertheless, as the equivalence also suggests, quantum-limited heterodyne detection does not deteriorate the 
sensitivity with respect to frequency-independent homodyne detection, except for the lower effective optical power, 
which can in principle be made as close as possible to the true optical power, as —> +00. For certain specially 
designed interferometers, such as the speed-meter interferometers i22j with Michelson [2^ or Sagnac J9,.30J topologies, 
the optimal homodyne angle is largely constant over a broad frequency band. These interferometers already exhibit 
broadband QND behavior with frequency-independent homodyne detection. In this situation, a heterodyne detection 
scheme (e.g., the Schnupp square- wave demodulation scheme), optimized for that particular quadrature, can be 
employed, e.g., for technical reasons, without compromising the sensitivity. 

V. CONCLUSIONS 

In this paper, we applied a quantum optical formalism to a heterodyne readout scheme for advanced GW interfer- 
ometers such as Advanced LIGO. Our results provide a foundation for the astrophysical optimization of Advanced 
LIGO interferometers and should be used to decide whether a homodyne or heterodyne readout scheme is more 
advantageous. 

One of the advantages of the heterodyne readout scheme (with the exception of balanced heterodyning), is that all 
output quadratures are available for measurement, providing a way of optimizing the sensitivity at each frequency. 
This result cannot be easily achieved in homodyne detection. Howe ver, as originally discovered by Gea-Banacloche 
and Leuchs |0 and by Schnupp ^3 1 analyzed by Niebauer et al. 0| , and Meers and Strain in the low-power 
limit, heterodyne detection leads to an additional noise term which is a direct and necessary consequence of the 
Heisenberg uncertainty principle. 

In the specific case of detuned RSE interferometers planned for Advanced LIGO, we derived the expressions for 
the total heterodyne noise spectral density [see Eqs. H22|l - H26() . Hill) ], assuming a pair of Schnupp sidebands with 
arbitrary amplitude ratios. In the balanced case the effect of the additional heterodyne noise is shown in Fig. 01 In 
the more practical very unbalanced 16] configuration, we compared the noise curve in the optimal heterodyne case, 
obtained by maximizing over the heterodyne phase at each sideband frequency, with some noise curves obtained when 
the homodyne readout scheme is used. The results are shown in Fig. [S] Neither the homodyne nor the heterodyne 
readout provides a noise spectral density that is the lowest for all frequencies. Moreover, the differences between 
the noise curves occur mainly in the frequency band 70-200 Hz where other sources of noise in Advanced LIGO will 
probably dominate, e.g., thermal noise (unless more sophisticated techniques are implemented |2^). So, before 
drawing any conclusion on which readout scheme is preferable, the comparison between them must take into account 
the other sources of noise present in Advanced LIGO and should be addressed with reference to specific astrophysical 
GW sources, such as neutron-star and/or (stellar mass) black-hole binaries, for which the GW spectrum is a power 
law with an upper cutoff ranging from ~ 200 Hz to several kHz, and also low-mass X-ray binaries which require 
narrowband configurations (detuned RSE) around 500 — 700 Hz. In this paper we have provided a framework in 
which these optimizations can be carried out. We shall report on the results of the optimization elsewhere 16]. 

From a more theoretical point of view, we worked out a frequency-independent quantum limit for the additional 
heterodyne noise [see Eq. H51|l ]. which made more explicit the following fact: lowering the additional heterodyne noise 
while simultaneously retaining the ability to measure more than one quadrature is incompatible in heterodyne detec- 
tion, which is inherently frequency-independent unless frequency-dependent squeezing techniques are implemented. 
In particular, this incompatibility seriously limits the extent to which conventional interferometers can beat the SQL 
using heterodyne readout scheme. Indeed, we show in Sec. IIV Bl that conventional interferometers with quantum 
limited heterodyne detection are equivalent to conventional interferometers with frequency-independent homodyne 
detection and lower optical power. However, for third-generation GW interferometers with speedmeter-type configura- 
tions [23, in, |2^ |33] , which are already QND interferometers under an appropriate frequency- independent homodyne 
detection, heterodyne readout schemes can in principle be employed without compromising their sensitivity. 
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